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Abstract 

The equations of motion of charged particles under the influence of short 
electromagnetic pulses are investigated. The subcycle regime is considered, 
and the delta function approximation is applied. The effects of the self force 
are also considered, and the threshold where radiation becomes important is 
discussed. A dimensionless parameter is defined that signals the onset of 
radiation reaction effects. 

Very short and even subcycle optical pulses have been gaining increasing 
attention in recent years. [Tj The theory of the interaction of short pulses with 
charged particles has been studied in [2|, [H] in one dimension. In three dimensions 
the subcycle problem becomes more complicated, and contentious jSJ-jHl; and it has 
also been studied in plasmas. [TU] Often, as the pulses decrease in their temporal 
span, the intensity rises correspondingly. In fact, intensities of 10 22 W cm -2 have 
been reached, and this number is expected to go even higher. JT] At such extreme 
conditions, the radiation reaction force should be examined, and we will examine the 
intensities and pulse durations where the onset of the self force becomes important. 

To begin, a nonrelativistic approximation is used to assess the use of a delta 
function to model a short pulse. It is shown to be in agreement with the exact 
solution of a Gaussian pulse, the approximation improving as the Gaussian pulse 
becomes smaller. It is then shown how the delta approximation works for the 
relativistic case, and finally, self force effects are considered. 
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For a short pulse the slowly varying envelope approximation fails, and for a 
subcycle pulse the whole notion of a wave train is derailed. Opposite the limit of a 
monochromatic wave sits a delta function, and in the following we examine the 
usefulness of this, other extreme, approximation. The value of this approximation 
lies in the following observation. As the pulse becomes ever smaller, the exact pulse 
shape may not be known exactly, and furthermore, sometimes the detailed motion 
during the fleeting moment of interaction is not of interest anyway, where the final 
velocity and energy are of more interest. 

The most drastic approximation to be made is the one dimensional model, 
meaning that the pulse, or waveform is of the form f(z — ct), i.e., the plane wave 
approximation (although not necessarily a wave). Real pulses are focused down in a 
variety of ways and, regardless of how they are made, the field is a function of all of 
the coordinates. This has been discussed extensively in the references given above. 
Nevertheless, there are times when the 1-D approximation is useful and sheds light 
on the more exact physics. This is one of those cases, where we are able to examine 
both the delta function approximation and self force effects with relative ease. 

For example, consider Gaussian pulse. Non-dimensionalized units are used 
where z/L — > z and ct/ L — > t, where L is taken to be of the order of the wavelength 
of visible light (for numerical calculations below L = 5, 000A), and the field is 
polarized in the x direction. 



where w is a dimensionless parameter fixing the width of the pulse. This must be 
accompanied by a magnetic field 



For nonrelativistic dynamics one may ignore the magnetic field and integrate 

F = ma = qE, for a particle with charge q and mass m. It is helpful to define the 

impulse by 



E = Ee- ({z - t)/w)2 x 



(1) 



B = Ee-^ z - t)/w)2 y. 



(2) 



I = Fdt = qEwL^/c, 



(3) 



J —oo 



With this, integrating the equation of motion yields 
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-=2( 1 + erf( — ») (4) 
This is plotted in Fig. 1 at z = for different values of w. 

As w gets small, we may think of the force on the particle as impulsive, and 
in fact, approximate the Gaussian by a delta function 1 E = N5 x where 
normalization N is fixed by assuming 

Ee-WM'dt = J N5 {^) dt (5) 

which gives N = \/ttE. This insures that the impulse is the same in each case. The 
velocity is easily found. Assuming that the initial velocity is zero gives 

» = («) 




Figure 1: mv/I vs t for w = 1,1/2,1/ '4, showing that the smaller w is, the more 
steplike is the response. 

This shows that replacing the Gaussian by a delta function reproduces the 
solution of the equation of motion in the limit of a narrow pulse. In the limit that 
w —>■ these results show that the Dirac approximation precisely coincides with the 
exact result given by the Gaussian approximation. In reality, short pulses have a 
complicated profile, but these results show the advantage of the delta function 
approach: It changes the problem from solving a differential equation to an 
algebraic equation. 

To consider the relativistic situation consider the equation of motion 



1 Since this is a function of z — t it satisfies the wave equation, which also follows from the fact 
that = 
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d -f = J-F-v.. (7) 
clt mc 

We consider an electromagnetic wave of the form 

E = Ef(z-t)x (8) 

and 

B = Ef(z-t)y. (9) 

This represents a plane wave of amplitude E, polarized in the x direction, described 
by any dimensionless function /. With this (J7J) gives, 

dv° 

^ =«/(«"- « 3 ) (11) 

^ = aV (13) 

where the dimensionless parameter a = eLE/(mc 2 ). 

As these equations stand, the delta function approximation fails. This is 
because, since the velocity is like a step function function, the integral of (fT3|). for 
example, is difficult to assess. 2 A better way to proceed is to note that (JTUJ) and ([T3]) 
imply, 

v° = 1 + (14) 

which leaves the pair 

^ = .V (15) 



2 For example, / 8(t)S(t) = 1/2. Since w 1 is not exactly a step function, this result is inapplicable 
as well. 
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These imply, 



(16) 



v°=l + (^) 2 /2. (17) 
Using the integral of (|14j) in the right hand side of (|16|) yields 

v 1 = a J f(-r)dT, (18) 

a fascinating result. It states that the x component of the four velocity is essentially 
equal to the nonrelativistic three velocity evaluated at the proper time r. With this, 
(pPfj) . and ()17|) . the relativistic solution is completely determined in terms of the 
nonrelativistic solution for an arbitray ID wave form f(z — t). 

For example, the delta function approximation may be used in (fTHj) . With 
the integral of (JHj), and letting / — > y/7i5, we find 

v 1 = a ^6(T) v 3 = Tia 2 /29(T) v° = 1 + vra 2 /2^(r). (19) 

It is easy to see that this agrees with the asymptotic form of the analytical solution, 

w 1 (r) = ^(l + erf(r)) l (20) 

and the other components are found by simple alegra. For example, for an intensity 
of 10 18 W cm -2 , which corresponds to a = —2.7, the exact analytical solutions are 
plotted in Fig. |21 The asymptotic region is seen to agree with the delta function 
approximation. 

Another useful observation from (J19)) is that for high intensity, v 3 ~ v °. In 
particular this is valid for a » 1. For example, at I = 10 19 W cm -2 the difference 
between v 3 and v° is less than one percent. 

One may also solve (|17|) with (J2(Jj) for the time t in terms of the proper time r, 



2 \ I 2 

e~ T \ „ , /2 / r— \ 2e~ T 



t = -a 2 7r ^rerf(r) 2 + 2 \r + -j=rj erf(r) + r - ^-erf (V2r) + \ , (21) 

which is used in Fig. |3] 
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Figure 2: The asymptotic regions agress with the delta function approximation, for 
V vs. t. 




Having an analytical expression for the velocity is useful for looking at the 
self force. The equation of motion with radiation reaction forces is 3 

= af^v, + b(v a + v v v v v a ) (22) 

GST 

where the dimensionless parameter b = ct /L, t = 2g 2 /(3mc 3 ) and the overdots 
imply differentiation with respect to the proper time. This is called the Lorentz 
Abraham Dirac equation, but the other form of the equation of motion is obtained 
by starting with but using 

- = a /- 1V (23) 
in the terms multiplied by b, which leads to, 



3 We assume t> cr u <T = 1. In the literature, some take v a v' T = — 1, which changes signs in the self 
force. 
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dif 

— = af^v, + abr»Vp + a 2 b(rf^% + .TV* V"), (24) 

which is called the Landau Lifshitz Rohrlich equation. The LAD equation fell under 
bad times to the the runaway solutions, preaceleration issues, or the apparent need 
to invoke non-zero size particles, ^2] while the LLR equations seems to be only an 
approximation. For a nice introduction to the issues, with historical notes, one may 
consult Rohrlich 13], which has many of the seminal references and discusses the 
distinction between "self forces" and "radiation reaction forces." More recent work 
considers the problem in various dimensions, the effect of the magnetic 

dipole,[in][IZI, connections to QED[18 and vacuum polarization. |2"U]. mass 
conversion |19j. and hydrogenic orbits. [2~T] 

Whenever b is small, which is true for all but the most extreme light, it is 
sensible to expand the solution in terms of this parameter, 

v a = v a + iv a + ... (25) 

where ov a is the solution without self forces, \v a is 0(b), and so on. Then, for the 
plane polarized field used above we have, 

^ = af 1 v 1 +g\ v) + O(b 2 ) (26) 
^ = af( lV ° - lV 3 ) + <?V) + Oib 2 ) (27) 

d ^ = af 1 v 1 + g\,v) + 0(b 2 ). (28) 

where g a represents the self force, and where qV g is known, given by (fTBj) . ()14|). and 
(fTTj) . These equations show that v° = v 3 + 0(b 2 ), a very useful result, which allows 
us to use v — 1 = v 3 in the 0(b) equations. With this, we have, calling = o^ 1 and 
dropping the subscripts, to 0{b) 

dv° 

r- = afv 1 + ab{<Pf - af<f/2) (29) 

and 
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dv 1 

— = af + ab(f - a<pf) (30) 

It is noteworthy that these equations are obtained by using either the LAD 
or the LLR form for the self force. We use 

-({z-t)/ W f e -(r/w) 2 

/=- >~ , (31) 

w w 

where the last part holds to 0(b 2 ). One can see that the self forces become 
important as r = c?b\pK '/w —>■ 1. We examine this, and the accuracy of the 
approximations, for specificic cases. For an intensity of 10 23 W cm -2 (a ~ —850) and 
w = 100, b ~ 0.5, the velocities are given in Figs. EJandEl It is seen explicitly, as we 
expect, the effect of the radiation reaction is to reduce the final velocity of the 
particle. If the intensity is increased by one order of magnitude, so does r, and the 
results diverge drastically, as expected. 4 

As a final example let us consider a pulse of soft x-rays, where we take 

-((z-t)/ W ) 2 

f = cos(Q(z-t)) (32) 

where the dimensionless Q determines the frequency. In this case the delta 
approximation is invalid, but the expansion in terms of b is still useful. To zero 
order in b we have, 

»" = 4 oe 4 ^ r L - — ) + erf u + — )) • (33) 

This of course is real, which is seen directly by writing the error function in terms of 

the imaginary error function. Before considering x-rays it is interesting to consider 

the case that Q — 1 = w at 10 23 W cm -2 . This is just below the radiation reaction 

"threshold," but it is interesting to see how strongly relativisitic the solution is. 

This is evident in Figs. Eland[7[ which show the four velocity as a function of proper 

time and the corresponding velocity (dz/dt) versus t. 

For soft x-rays Q ~ 100 (50A radiation) and we take w — 1/10 (~ 200as 

pulse). The results are plotted in Figs. Hand El for / = 5 x 10 23 W cm -2 . 

4 As a partial check on the numerical work, which was accomplished using Mathematica, v a v a 
was plotted, a result not used explicitly in the calculations. The value was always within a tenth of 
one percent of unity. 
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Figure 4: vq plotted against lOOt with no ra- Figure 5: V\ 

diation reaction, as a solution to (|29~j) and 
f)3(jp . and a numerical solution to the exact 
equations |2H 




500 1000 1500 2000 



Figure 6: v z plotted against the proper time Figure 7: dz/dt vs. t, corrsponding to Fig. 
for I = 10 23 W cm" 2 . El 




Figure 8: vq plotted against lOt with no ra- Figure 9: V\ 

diation reaction, solution to (j2HI) and 

(pffij) . and a numerical solution to the exact 
equations EH 
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In order to assess the vailidty of the approximations we may start with the 
equations, for any pulse of form considered above, 



^! = a (f + bjy + a 2 bf(v° + v 3 - Z 2 v°) (34) 

GST 



^ = a(f + bf)Z-ba 2 f 2 Z 2 v' (35) 



^! = a (f + + 6a 2 /V + 1, 3 - ZV) (36) 
ar 



where Z = i>° — v 3 . From and (}3l)j) we have 



r 2 = 1 + 26a : 



/ frfr. (37) 



This shows that for b = 0, Z = 0, which we found above. Since the integral is 
bounded, Z « 1 for small 6. We can use this to investigate a self consistent solution 
to this equation by using the Z — 1 (which implies t — z = r) in / to find 

Z~ 2 = 1 + ba 2 w^ (l + ed(V2r/w)) . (38) 



This shows that the approximation Z = 1 is valid as long as ba « 1- 

In summary, it has been shown that the delta function approximation is 
applicable for both the relativistic and non-relativistic case for any incident traveling 
wave of the form f{z — t). For the relativistic case, it is shown how an exact 
solution may be found, and these results are used in the investigation of self forces. 
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